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ABSTRACT 


This work deals with characterization of deep level 
impurities, particularly of Gold in silicon. Isothermal 
relaxation data, using space charge spectroscopic principles, 
has been cdjtained, using set ups to measure lew level capa- 
citance end current signals from the device. Both single 
shot and synchronous detection methods of measurement and 
analysis haye been used. The problem of det, of multi^ 
exponential transients, which has been eluding an accurate 
and reliable solution, has been treated using an appropriate 
signal processing algorithm. In the present work the 
acquisition and analysis functions, v\hich are performed 
simultaneously in the classical deep level transient spectro- 
scopy (DLTS) developed by Lang, have been separated. This 
separation has been achieved using transient recording and 
analysis on a signal analyser, whoso large memory as well as 
facilities for development of various analysis and signal 
processing routines have been found of immense use in analysi 
The results have been obtained on a commercial p'^'N Jn 
diode 1N4148, for which the other relevant data was readily 
available. Important parameters such as activation energy 
and capture cross-sections have been determined. 
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CHAPTER 1 


INTRODUCTION 


Deep levels in semiconductors play an important 
role in device perfornence. They determine crucial device 
parameters as mean lifetime, activation energy, recombi- 
nation efficiency etc. Hence it is essential to accurately 
characterize deep levels. This, however, is not a simple 
process, essentially because deep levels can get created 
not only when a known impurity is added but also they get 
inadvertently created due to lattice defects. The presence 
of one can affect the performance of the other deep 
levels as a recombination center. For instance Gold is 
extensively used in semiconductor industry for fast swit- 
ching operations- and lot of data on its parameters has been 
reported in literature. The reported values, however, 
vary by orders of magnitude, Vlu and Peaker [l^] have 
reported a survey of these values over the years. Table 1.1 
just lists the electron cross-section values of Gold 
acceptor in silicon. This clearly suggests a degree of 
unpredictability related to values of Gold parameters. 

This unpredictability is partly j-iustified by Bullis [9], 
wherein these parameters are found dependent upon several 
factors like the characteristics of initial neterial. 
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Table l*l: Electron cross-se 
silicon. 


Author 

Year 

^nCcm^l ^ . 

Bemski and 

1958 

5x10“^^ 

Struthers 

Davies 

1959 

ZxlO”^^ 

Fairfield and 

Gokhale 

1965 

8. 4x10“^’^ 

Senechal and 

1968 

Z.SxlO"*^^ 

Bas inski 

Sah et.ai. ' 

1970 

2X10“^^ 

Zohta 

1972 

4,6x10“^^ 

Kassing and 

1974 

8,3x10“^^ 

Lenz 


1.3x10“^^ 

pals 

1974 

Nagasawa and 

1975 

2,2x10“^^ 

Schultz 


—14 

7.3x10 

Kasaing et.ai* 

1975 

Barbolla 

1976 

8. 3x10“^*^ 

et.ai', 

Dudeefc and 

1977 

1.5x10“^^ 

Kas'sing 

Brotherton and 

1978 

9,10"^'^ 

Bicknall 

Nassibian and 

1978 

1.10*'^^ 

Farone 

Lang et.ai 

1980 

6,9x10'“^’^ 

Lang et.ai. 

1980 

1,7x10’*“^ 

Wu and 

1981 

8.5x10“^*^ 


Peaker 


;tion of Gold acceptor in 

Te mp(K) Material Technique 

200-450 Nj=N^ Minority 
^ Lifetime 

77 . CZ NT,»Nn Majority 

^ Lifetime 


300 CZ Nt^O. 1 P,C. Decay 


250 

CZ N^Nq 

June. Cap, 

80 

N^O.a^P 

Maj. Flux 
junc, cap. 

300 

FZ N^Nq 

June, Cap, 

300 n+i n'^N^>Njj 

SCLC 

77-220 

CZ Nt=0.2 

■"d 

Puis e/Junc . 
cap. 

300 

Implanted Puise/Junc, 
Cap. 

172-316 

n"^! n‘^N^>Nj^ XLC 

80-205 

CZ Nt^O.3 

Nd 

Pulse/junc. 

Cap. 

300 n'*’i n'^N^>NQ 

SCLC 


77-280 CZ N^. IN Q Puise/junc .C£ 

235-265 Minority/1 if el 

295 EPI N^<0,1Nq Pulse/Junc.Cai 
115-424 CZ N^O - IN Pulse/ June, cj 
80-200 CZ N^ O.lNp Pulse/Junc.ci 
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behaviour of one type of deep level in presence of another 

deep level or even the shallow levels. Lang [ 10 ]bas conclude 
that in place of a well defined unique level, there -are family o 

closely related but different gold based complexes. 

This work concentrates on Gold in silicon. The deep 
level inpurity characterization can be done using either 
isothermal or nonisotherraal relaxation. Here the former 
method is chosen and in particular, junction spectroscopic 
methods have been used. These methods have the advantage 
that one can distinguish between electron and hole traps 
and further one can get directly, emission rates and related 
data. Both single shot and synchronous methods of analyses 
have been used. They are detailed in Chapter 2. Chapter ■ 

2 also includes the theory behind current and capacitance 
transients as tools to deep level characterization using 
isothermal junction transient methods. Lang's method of 
Deep Level Transient spectroscopy (DLTS) is also explained 
in Chapter 2, 

Chapter 3 discusses the analysis procedures for 
current and capacitance transients from the sample. It is 
found that this data is quite unamenable to accurate 
characterization because of two reasons, one that the 
transient has a steady state level on which it sits and two, 
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that instead of having single exponent, it can have a sum 

of exponents (ox multiexponent). Physically multiexponents 

can mean that more than one trap are involved in recombination 

process. The data in that case can thought of having 
n . 

form like Z A- e““’i , The problem is to accurately 

i=l ^ 

find values of n and Apparently it looks a trivial 

problem, which it is not [ll]* According *to Lanczos [11} , 
the practical solution of these problems is beset with 
numerical difficulties. For an accurate deteraination , large 
data accuracy is needed. If it is not so, one can get 
totally unsatisfactory results. In no v\ay, any of the least 
Square treatments can help retrieve the situation for insu- 
fficiently accurate initial data. Chapter 3 goes into details 
of these problems, Ji/kiltiexponents have been extracted using 
an appropriate signal processing algorithm due to Gardner 
et.ai, [l2j. Chapter three contains the details of this 
method and the reasons of its success. In the same chapter, 
several other signal analysing algorithms axe described. 

They deal with the determination of parameters in a single 
exponential decay from the samples, in the presence of a 
constant d.c, shift, > 

The : particular method of . 

DLTS suffers from the fact that analysis is done using 
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hardwired instruments and each point in Arrhenius plot 
requires one thermal scan. Obviously this is very time 
consuming and hence this work delineates it into two inde- 
pendent functions of acquisition and analysis. The details 
of Box-car programs in Appendix I, help in underlining the 
advantage of using one thermal scan and getting all what 
one could get from conventional IXTS. 

Chapter 4 contains the de-fcails of measurement system 
used for characterizing the deep levels. The design of 
capacitance and current transient measurement modules is 
described in the same chapter. In both the cases ^ the signal 
is of low level in nature. The capacitance transient is, 
in particular, monitored by high frequency sneil signals. 

Chapter 5 contains the experimental data taken on a 
commercial p'^N junction diode 1N4148. Using the theory and 
analyses of second and third chapters, deep level parameters 
like activation energy and capture cross-section have been 
determined. Discussions of the results follow in the same 

chapter. 

Chapter 6 contains the conclusion and suggestions 
for further improvements , 



Appendix I and II contain program listings for 
analysis procedures described in Chapter 3. 

As lot of data has been obtained from the signal 
analyser, a typical plot is explained briefly in 
Appendix III, 


CHAPTER 2 


ISOTHERMAL RELAXATION STUDY RDR DEEP LEVEL IMPURITIES 

2.1 INTRODUCTION: 

Current transport in semiconductors is determined 
by electrais and holes in conduction and valence bands 
respectively. The number of electrons in the conduction 
band or holes in the valence band is dependent upon the' 
doping and the generation-recombination process in the 
semiconductor. At any instant these carriers can be 
generated either thermally or by any other, external per- 
turbation and they recombine with a mean life time. The 
process of recombination, particularly, is strongly 
dependent upon the impurities in the semiconductor. Hence 
efforts have been made to measure the parameters of the 
impurities (e.g, life time etc.) accurately and relate 
them to various device parameters like conductivity etc. 
The process of recombination can be band-to-band or 
band- to- impurity- to band. However theoretical values 
assuming only band to band transitions have been orders of 
magnitude different from experimental results in silicon 
and Germanium, Shockley et.al, [6] in their theory 
called SRH statistics explain the discrepancy in these 
’ indirect* materials by assuming the existence of one or 
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more deep levels within the band gap. These levels act 
as recombination centres. Tliough this idea of the existence 
of deep levels in SR^ statistics was a theoretical 
construct, subsequent experimentation has confirmed their 
presence. Thus an accurate characterization of these 
deep levels has helped solve two problems. One that the 
above said discrepancy in values of mean lifetime has 
been explained and secondly one can tailor device chara- 
cteristics according to one's needs. The exact effect 
of deep levels on recombination can be known from SRH 
statistics. Here only those aspects of SRH statistics 
will be discussed which are pertinent to our characteriza- 
tion of deep levels. 

Consider a deep level, acting as a recombination 
center within the band gap as shown in Fig. 2.1. The deep 
level recombination center or trap can essentially; 

1, capture electrons from conduction band r^ 

2. emist electrons to conduction band rj^ 

3, capture holes from valence band r^ 

4. emist holes to valence band r^^. 

Let <V> be the mean thermal velocity 

o^jCTp capture cross-sections for electrons 
and holes 

ttie trap. concentration 
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e ,e^ the emission rates of electrons and holes 
n' p 


and 


f = 


l+e 


'(“E be the probability that 


a given trap level is occupied by an electron. 

The rates of the above four processes 3:g,r|^,r^ and r^ is 
given by 


1. 

^a = 

<V>0^n*N^(l-f) 

(2.1) 

2. 

^b = 


(2.2) 

3. 


<V>ap.p.N^, f 

(2,3) 

4. 

^d = 

ep.N^(l-f) 

(2.4) 


under thermal equilibrium 
Ta = and r^ * r^ 

From Xq = we have from eqn, (2,1) and (2.2) 


<V> aj^*n.N^(l-f) N^.f 

on solving eqn. (2.5) one gets fox e^ 

(E+-.E )/kT 

e„ = <V> a„ N, e 
similarly from r^ = r^^ one gets 


<V> o^ n e 


(E -E.)/kT 


v 


(2.5) 


( 2 . 6 ) 


(2.7) 
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From eqns, (2.6) or (2,7), if one plots natural logarithm 
of e or e against l/T, the resulting strai^t line gives 

1 1 JH 

the activation energy as slope and ordinate intercept 
provides capture cross-section. Such a plot called Arrhenius 
plot is used to characterize the deep level recombination 
centre. Hence methods, both direct and indirect, have been 
devised to get values of e or e . In this analysis direct 
method using junction spectroscopic principles has been 
utilized. It has been found that the capacitance or 
current due to the deep level filling or emptying in the 
space-charge region is in the form of a transient and is 
related to emission rates of carriers and the . trap 
concentration. Section 2,2 gives in detail the origin of 
the capacitance and current transients. The analysis of 
these transients can be done in several ways and are 
detailed in Section 2,3* 

2,2 CAPACITANCE AND CURRENT TRANSIENTS AS PROBES TO 
RELAXATION STUDY: 

Let n^ and Pj be the concentration of electrons and 
holes at the trap. From eqn, (2,l) and (2.2), the rate of 
change of electron concentration (^)-j in conduction band 
due to trap can be found as 

= ®n "t - <=n " Pt 




(2.S) 
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where = <V> 

Similarly from eqns, (2*3) and (2.4), for rate of change 
of hole concentration in valence band due to trap, 

can be written as: 

From eqns, (2.1), (2.2), (2,3) and (2,4), the rate of 

dnj 

change of trap electron concentration (-g^) i® 


("dt ) = ^ 

= W + 

+ {dp + V>*^TT 


c 

p 


)n 


T 


( 2 . 10 ) 


where = n^ + Pj = trap concentration. Solution of 
eqn. (2.10) yields within depletion region for n=p=o 


n^(t) = n^(~) + [n^(o)-nj(«’)]e 




( 2 . 11 ) 


where, 

nj(a») is the final value and is given by 

n^^oo) = ) 

n p 


( 2 . 12 ) 
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Hr: 

and nj(o) = ^ f O ' e "" ^ (2.13) 

X 

= (P% Jn) 

« 0 (W^P Jn) 

where Ng is the equilibrium concentration, From eqn, (2.11) 
one sees that trapped electron concentration decreases 
exponentially with time and hence within deplection region, 
the corresponding capacitance and current, too, will 
change according to trapped electron concentration, 

A typical experiment to measure this capacitance 
or current transient involves application of a suitable 
voltage ■ . - observing 

the sample response. The response, however, is strongly 
dependent upon temperature. Therefore one does the 
experiment either at a constant temperature and hence iso- 
thermal measurements, or varies temperature according to 
a program and therefore nonisothermal measurements. This 
work uses isothermal conditions and the applied voltage 
is in the form of a pulse-bias. There are two' types of 
pulses one can use. In one there is no carrier injection 
process in the sample and in another there is carrier 
injection. The former is called a majority carrier pulse 
while the latter is called minority carrier pulse. Bach of 
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the pulses can be used to characterize majority and 
minority carrier traps respectively [5], The effect of 
a pulse bias on filling and emptying process in P'^N 
junction diode is shown in Fig, 2,2. The pulse used is 
a majority carrier pulse. In subsequent sections, both 
capacitance and current transient signal expressions will 
be derived for a majority carrier pulse in a p'^N Jn. 


2,2*2 Theory of Capacitance Transients: 


The transition capacitance in a junction diode is 
given ^ ® ^ * From Poisson’s equation and its 

solution we get 



(2.14) 


where W is the width of depletion region, Nj the con centra - 
tion of ionized impurities in it, £ the permittivity within 
it and V. the junction potential. 

From eqn, (2,14) and , one gets 



(2.15) 


For an acceptor level in a P'^N junction, Nj- 
hence eqn. (2.15) becomes 


Nq - n^ and 


“2?r 


(Nj-j “ n^) 


( 2 . 16 ) 
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Fig. 2.2: 1) Majority carrier pule.® (2). (4) : Energy band diagr 
Showing depletion region t with shaded lines) as well 
emissim ©nd trap ^occypation during the 
_ ^ c of the nejority carrier pulse 
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(2,17) 


(2.18) 


(2.19) 

carrier pulse 

gives value of 'Z^ directly. Using eqn. (2.6), one can 
obtain values of activitation energy and capture cross- 
section. The initial value of transient has informa ti<ai 
regarding the trap concentration. 

2.2.3 Theory of Current Transients: 

The pulse sequence applied in this case, too, is 
as in Fig. 2 , 2 (l).As the voltage switches from 0 to 
three components are obtained in the current transient 
signal. 


or 

^T 

Where, 


c, = [ 1 - Aclli ] 

m 


m 


qN, 


and 


AG(t) 


m 


n- 


Using eqn. (2.13) a and for acceptor level 
get for (2,18) 


C(t) _ 


m 


2N, 


N 


IT -Vtn 


= 3t; ® 


Thus capacitance transient from a majority 
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1. A fast component in the foxm of a spike of very- 
short duration is produced due to svs/eep out of majority- 
carriers in the part of depletiOTi region, which gets 
affected due to pulse bias sequence, 

2, The slow current transient due to trap occupation 

has two forms. One consisting of conduction due to 
electrons and holes and other due to displacement currents, 
because of induced charge in x < o and x > w region, due 
to • The components of slow transient are shown in 

Figure 2,3, 


The Gcaiduction current due to electrons 


w 


A / q (^)x dx 


n 


A q e^n^(t) w( t) 


( 2 . 20 ) 


Where A and vj axe area and width of junction respectivel-y 
and n=o while using eq. (2.8), 

3, The displacement current is given by 

w dn~ 

i.(t) - A / q( ^ ) (x/w)dx 


dn. 


qwA /'"‘T 

“ 2 '"W 


) 


( 2 . 21 ) 



V 

d= t= to, 


FF 


jk. 




e- 


Fig* 2*3: Energy band diagram during the svvitching transients 





Combining eqn, (2.20) and (2.2l) and ignoring conduction 
current due to holes, which will be small in n-type 
semiconductor: 


i =» ®n (2.22) 

and for case of P'^N Junction (2-22) becomes 

1 = 2|A (2.23) 

For an acceptor level (2.23) becomes 

i « e (2.24) 

using the expression for V/ in (2.15) and (2.24) becomes 

i _ ^ __ %T j-l/2 

2 Mjj 


(2.25) 


where. 


W 


[ 


3 1/2 


q Nq 

If << Np, (2.25) becomes: 


qAejjNlT 

2 


W e 

OO 




Hence the current transient is exponential when N 
otherwise it is ncaiexponential. 


(2.26) 

TT ^D* 


The timing diagram fox both current and capacitance 
transient in case of majority carrier pulse is shown in 
Fig. 2.4, 

2,3 TECWiqUES USED IN ANALYSIS OF TRANSIENTS; 

This Section briefly discusses the various techniques 
used to get trap data. They are essentially of two types 
called single shot methods and synchronous detectim methods. 
In both the aim is to get Arrhenius plot, 

2.3.1 Single Shot Methods: 

They were historically the first methods to be used 
for analysis and are accurate but axe very painstaking. 

Prime examples of this in our case are conventional iso- 
thermal capacitance transient study and isothermal current 
transient study [4], Herein ■Uie time constants of transients 
are found and each point in Arrhen in' s plot is obtained 
thereby, 

2.3.2 Concept of ELTS: 

Deep Level Transient spectroscopy is a synchronous 
detection method. It has been used very extensively in 
getting trap parameters. The major feature of this method 
is that it is fast, sensitive and spectroscopic in nature 
[5], A wide variety of traps can be detected by observing 
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the DLTS spectrum. The essential feature of IXTS is that 
one can set a rate window in a Box-car integrator as 
depicted in Fig, 2,5. The difference signal C(t^)-C(t2) 
or i(t2^)-i(t2) is plotted as a function of temperature. 

The sign of the trap tells whether it is an electron or a 
hole trap and temperature corresponding to maximum gives 
one point on the Arrhenius plot. The value of ®p) 

needed can be related to difference (t^-t2) through following; 

Normalised DLTS signal is given by 

= [C(t^) - C(t2)]/-^^C{o) 

where C(o) is the capacitance change at t = o. 

For exponential transients 

“t-i/'c. 

S(t) = [ e 1 - e ^ ] (2.27) 

Where is the time constant of the transient. 

Maximising eqn, (2.27) with respect to by differe- 
ntiating and equating it to zero gives 

= (Vtp [in [Vtg]] (2.28) 

Thus one needs simply to scan the temperature and observe 
the K-TS spectra for each scan. The sign of peaks and 
position of peaks need to be observed and from this and 
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values of and using eqn. (2*28), one can get an 
Arrhenins plot. Some remarks are in order: 

1, This method is applicable to exponential transients 
only as evident from eqn, (2,27). Any non-adherence to 
this criterion will give wrong values for activation energy, 
capcture cross-section and related parameters [7], It has 
been further -found that choice of t^ and t^ can give widely 
different Arrheni s plots in such cases, 

2 , Even for a single exponential, the conventional 
DLTS uses inefficiently the data, it gets in each scan. 

Hence modifications have been suggested [8] where the rate 
windowing concept is done in softv/are after recording all 
the transients in a single scan, 

3, DLTS subtracts the offset due to leakage current 
or steady state capaci'tancG and thus overcomes a serious 
problem affecting the single shot analyses, 

4, Finally DLTS is at best a very good survey technique 
and can in no way supplant careful and accurate 
determination possible through single shot analysis. 



CHAPTER 3 


TRAI-^SIB'IT DATA AI'JALYSIS 


3.1 INTRODICTION; 

The isothermal relaxation data is characterized 
by two major factors; 

a) That the transients sit on a steady state level 

A Iq oi -^^ 0 ^ leading to the form + B. 

b) That in the presence of multiple traps, the data 

-a t 

has multiexponential decay of the form Z A-e i 

i=l ^ 

At the outset, one might feel that determination of 
the values of A^ and a in the above is a trivial problem. 
However, it is not so [ll]. Even for a single exponential 
case in (a) above to get an accurate and reliable value for 
A, a and B is a difficult exercise. In this chapter liiese 
aspects of the problems have been studied, 

3.2 ANALYSIS FOR SINGLE B<PONENTIAL: 

Initially the problem of determination of decay 
constant in single-exponential transient seems very easy. 
One can "Wiink of simply finding the ratio of the transient 
curve at two values and computing the decay constant. Or 
one could fit the logarithm of the transient curve to a 



26 


straight line using linear regression algorithm, thereby 
getting both A and a. However the offset B frustrates any 
Such simple implementation. Hence the aim is to somehow 
remove this B and use the conventional methods thereafter 
or devise a suitable algorithm v/hich takes B into account. 
The problem is further exacerbated by certain limitations of 
the signal analyser, detailed later on in the chapter. In 
any case, one has to be careful. Several methods have been 
suggested to remove this B accurately, particularly in case 
of isothermal relaxation 'tvo ivays were found to determine 
the value of B: 

a) One uses the physical origin of this constant 

offsets In either current or capacitance transients, AIq or 
^ emanate from the leakage current or steady state 
junction capacitance during the time of reverse bias in both 
minority carrier and majority carrier pulses. Hence for a 
given device, one can reverse bias it to the same conditions 
as that of pulses and obtain values of AI and Howeve 

the major problem in this lies in obtaining identical condi- 
tions repeatedly for each tem.perature. Hence it is not a 
very attractive method, 

b) The other method to determine B, is to take 
average of a number of points on the tail of the transient 
and then subtract this average value from the whole curve. 
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However, scatter due to noise in the signal precludes 
dependence on only a few points. Hence c«ie needs to take a 
large enough number of points on the tial portion so that 
randomness of noise cancels and helps in accurate deter- 
mination of B, Such requir einent, particularly, in case of 
low temperature slew transients is not readily possible. 

One can then think of logical extension of the basic 
methods. One such extension is: 

3.2,1 Ratio Method: 

Herein we take three points on the transient and 
solve the resulting transcendental equation. Its brief desc 
ription is as follows: 

Three points on the transient curve + B are: 


-ati 

= Ae + B ( 3,1) 

= Ae ^ + B (3.2) 

-at- 

73 = Ae + B {3.3) 

From ecpis, (3.1), (3,2) and (3.3), we get 

-at- -at, 

72 - Yi = ■^(e ^ - e ^ ) (3.4) 

-at- -at- 

73 - 72 = A(e - e ^ ) (3.5) 
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Using eqn» (3*4). and eqn, (3,5), one gets 


^ 3-^2 


“at« -at, 
e ^ - e 



(3.6) 


Values of Yi»V2*'^3* ^1*^2 ^3 known from the transient 

data and one has to solve this transcendental equation. Using 
Newton-Raphson method, one can solve this iteratively. 


Normally one should get a reliable and satisfactory result. 
However numerical snags come up, \vhen evaluating it. 

This is because even slight variations in the values of t^, 

^2 ^3 widely varying values of a^s. Actually this 

is symptomatic of all problons dealing with exponentials. 

This method, therefore,, demands very accurate measurements 
of tj^, t2 and Table 3,1 emphasizes the above ccaitention. 

An alternative method is, hence thought of: 


3,3,2 Area iviethod. 

Here transient with a positive d,c, shift is assumed, 
though for negative offsets too, the same procedure holds, 

A (current) transient fr<»i the relaxation process with a 
positive offset is shown in Fig, 3,1. 

The three equations in this case are: 


= A + B 


(3.7) 
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Table 3,1: Results from ratio method analysis 


^0 
in ps 

^1 

in ps 

‘‘^2 
in ps 

. ^3 
in pS 

A 

a 

B 

0 

0 

0 

0 

.0^285603 

13610.6 

0.226716 

+2p 

X 

X 

X 

0,^267275 

13621.6 

0.22678 

+2p 

X 

X 

X 

0.233989 

13631.7 

0.226838 

-2^ 

X 

X 

X 

0.305109 

13592.2 

0.226608 

-ep 

X 

X 

X 

0.348874 

13633.9 

0.226851 

0 

+2p . 

0 

0 

0.281506 

12557.3 

0.220042 

0 

+4p 

0 

0 

0.279169 

10867.5 

0.206545 

0 

+6p 

0 

0 

0.312894 

6533.76 

0.139164 

0 

-2p 

0 

0 

0.287909 

14063 

0.229265 

0 

-4p 

0 

0 

0.295239 

15257 

0.235239 

0 

-6p 

0 

0 

0.292871 

14884 

0.233479 

0 

-lOp 

0 

0 

0.300032 

15863 

0.237912 

0 

0 

+2p 

0 

0^291945 

16529.6 

0.242591 

0 

0 


0 

0.290297 

16061.2 

0.240467 

0 

0 

+ 6p 

0 

0.287189 

14965 

0.234946 

0 

0 

+iop 

0 

0.292163 

10658.4 

0.200881 

0 

0 

-4p 

0 

0.297631 

17829.1 

0.24783 ‘ 

1 

0 

0 

-6p 

0 

0,295431 

17367.7 

0,246071 

0 

0 

0 

+4p 

0.285849 

13564 . 

0.226259 

0 

0 

0 

+ 6p 

0.290505 

12782.7 

0.218111 

0 

0 

0 

+iop 

0.28221 

14338,3 

0.233456 

0 

0 

0 

-4p 

0.291403 

126505 

0.216632 

0 

0 

0 

-6p 

C. 293571 

12349.9 

0.213153 

0 

0 

0 

■■15p 

0.349857 

8266.78 

0.1405 

0 

+2p 

4-4p 

+iop 

0.287851 

15939.9 

0.240383 

0 

-icp 

-I5p 

+10p 

0.299729 

17945.7 

0.246896 

+3p 

. +5p 

-7p 

+2+ 

0.259482 

16174,4 

0.239715 

-2p 

.■*«5p 

4-4p 

+9p 

0.319644 

16026.8 

0.238732 

tl = 

43,^^r 

pS, 't^s: 

ps 


= d 


p = 2.4414 ^s, X =: Corresponding change, \f'^2 "■ 

in t^,t2 etc. 
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-at. 


= Ae 


and 


+ B 


-.at. 


K = / (Ae + B) dt 


(3.8) 

(3.9) 


eqn. (3,9) on evaluation, gives 

-at, 

Ka = A[l-e 3 + Bat^ 


(3.10) 


Us 


ing B = Xj^-fA from eqn . (3,7) in eqn. (3.8) we get for A, 


- ^2 


-at, 

( l- e ^ ) 


(3.11) 


Substituting eqn. (3.11) in eqn, (3.10), we get: 
(X^-X 2 )[l-at^-e"*“^ ] (l-e ) (K-X^t^)a 

(3.12) 

The above transcendental equation can be solved by observing 

the zero cross-over, for monotonically increasing a. 

B 

The values of ^ m obbS • obtained from this method 
have been found to be accurate and less susceptible to 
initial data variations than the Ratio method. Table 3,2 
illustrates the above point. It further shows that the 

variation in a value with different choices of t ( t=o) is 

, o 

not acceptable. Unfortunately an accurate determination of 
is not easily possible for current transients. This is 
precisely due to two factors: 
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Table 3;2 

Results ficxn Area Method 


"o’ 

■ i 

a 

B 

0 

O'- 

12825 

0.216796 

+2p 

0 

9775 

0.216796 

+6p 

0 

9675 

0.216796 

-2^ 

0 ■' 

13825 

0.216796 

-4p 

0 

. 9525 

0.216796 

0 

+2P 

13550 

0,217285 

0 

-2p 

14550 

0.217285 

0 

-4p 

14225 

0.218261 

0 

-6p • 

16025 

0.219726 

+2p 

.2p 

10100 

0.219726 

-3p 

-4p 

10000 

0.218261 

-P . 

+P 

9500 

0.217285 

>• 

^0 

0, tjj^ s 2,13378 ms, 

p « 2.4414 ^s 

Jill M .HI. I...■..l ■ r .L|.. . nr 


^0*^1 ^ ^ % re^ectivelyi 
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a) 100 KHz bandwidth of the signal analyser and 
hence existence of reasonably strong Gibb’s effect.- 

b) A large current component due t^ sweep-out of 
majority carriers in the sample, saturates the current 
transient measurement set-up. 

The cumulative effect of both is that accurate ■^o 
determination is not possible, Hov^ever the redeeming 
feature of the area method, as from Table 3.2, is that the 
variations in B, due to fluctuations in choice of t^ instant 
are very small. Therefore our analysis for single expo- 
nential transients with offset, starts with Area method, 
where frcHii the value of B is determined. This B is 
subtracted from the total curve and conventional methods are 
applied to get the values of A and a accurately. 

The above analysis proceeds on the assumption that 

a single exponential decay is present, Vs/hat if it is not 

so, as it is in the case of multiple traps where we may force 

n „ + 

unknowingly a multiexponential of the form S A. e i 

i=l ^ 

to a single exponential. According to Lanczos [ll], 
satisfaction of equivalence of the measured results to the 
computed ones using either one or more decay constants is 
no indicator of the actual values. For example the data 
in Table 3,3 for a double exponential signal can be, within 
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Table 3^3i Illustration of numerical equivalence 

of two different exponential expressions. 


t 

in 

microseconds 

Actual double 
exponential 

Least Square single 
exponoitial 

2,98375e”^°^^**^^^ 

0 

2,99988 

2.98372 

2,4414 

2.93957 

2,92562 

4,8828 

2,88049 

2,86863 

7,32418 

2.82275 

2.81276 

,9.76558 

2,76623 

2,75797 

12,2069 

.2,71094 

2,70426 

17,0897 

2.60388 

2,59994 

21,9725 

2.50134 

2,49964 

24,4139 

2,45178 

2,45096 

26,8553 

2,4032 

2,40322 

3^,1795 

2.26379 

2.26551 

39,0623 

2.17579 

2,17811 

43,9451 

2,09155 

2,09412 

51,2693 

1.97168 

1,97416 

56,152 

1,896 

1,89804 

63,4762 

1,78^5 

1,78937 

70,8004 

1,68701 

1,68697 

80,5659 

1,56201 

1,55961 

85,4428 

1,5034 

1,49966 

90,3315 

1,44702 

1,44109 

95,2144 

1,39307 

1,3855 

100.097 

1.34131 

1.33205 

Least mean s 

qua re error for 300 

points = .151819x10"'^ 


35 


the accuracy of measurements, fitted to a single exponential. 
Comparing the two data lists one finds that error indeed is 
very small. In order to be more confident in case of single 
exponential analysis as well as to ascertain the existence 

of multiple traps, it is mandatory to knew values of n, A. 

^ - t ^ 

and in E hi e i , This is dealt in next section, 

^ i=l ^ 

3.3 ANALYSIS FOR iCTLTIEXPONH'JTIAL SIGNALS; 

Several methods are knovm for this important analysis 
According to Lanezos [ll], it is not the limitation of the 
methods of analysis, viiich gives us unsatisfactory results 
but that the accuracy of initial data is not good enough 
to take care of extreme sensitivity of exponentials to 
scatter in initial data. Hence a method due to Gardner 
et.al, [12] is described. In this method the end result is 
in the form of peaks corresponding to each of the exponents 
in the multiexponential signal. Tliis method is subdivided 
into three main subsections viz, theory, implementation ? 
results and application to .relaxation data, 

3,3,1 Theory: 

n -a-t 

VVe are given the function f ( t) = Z A. e and 

' 1=1 ^ 

we wish to determine n, A and The above Dirichiet 

series is converted to Stieltjes integral; 
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f(t) = E A. e = /“ dh(a) (3.13) 

i=l ^ o 

where h(a) is a sum of step functions, Eqn. (3,13) can be 
expressed in terms of Delta functions as: 

fit) = / e“^^ g(a) da (3,14) 

o 

The quantity g(a) contains information about the various 
values. The aim is hence to solve for g(a). The 
solution of this integral equation is done by transforming 
the variables a and t as a = e“^ and t = e^. 

Then eqn, (3,14) becomes 

f(e^) = / exp[e“^^"'^^]g[e“'^]e~^dy (3.15) 

comparing eqn, (3,14) and (3,15), one finds 

g(a)da = ag(e--^)dy 
or 

(^^■^ )da = g(e”^)dy (3,16) 

Hence a plot of g(e‘‘^) with respect to y is equivalent to 
plot of g(a) with respect to a. Therefore we need to find 
ways to get g(e“^) from eqn, (3,17), We multiply ecp, (3.15) 
by t=e^, then eqn, (3,15) becomes: 
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f( = / exp[“e^“^3 g(e""^)dy (3.18) 

-^•OO 

Vie define Fourier transform F(ji) as 

F(p.) = /” e^f(e'')e^^’''' dx (3.19) 

Hence eqn. (3.18) becomes using eqn. (3,19), 

?(}!) = f S exp[~e^“'^3 e^“'^g(e“'^)dy exp(ipx)dx 

^ .-iCO —CO 

(3.20) 

Putting s = x-y and x = s+y in eqn. (3,20) one gets; 

F(}i) =-^~ / / exp[-e®] e® g(e‘'^)dy 

* . — oo. — oo 

exp( ip,(s+y))ds (3.2l) 

ox eqn, (3,2l) becomes 

F(p) 9(e“^) exp(ipy dy / exp(-e )e 

exp(ivis)ds (3,22) 

In eqn. (3,22), one finds that first integral can be called 
Fourier Transform G(p,) of g(e“^), while the second integral 
needs to be found separately, Hov-'ever the second integral 
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can be related to Euler integral of complex Gamma function, 
which is given by: 

P{z) = / e~'^ t^~^ for > o (3,23) 

wCO 

Substituting t = e® in eqn, (3,23) one gets 

r'(z) = / exp(-e®) e® ds (3,24) 

••OO 

If z - 1+ip,, then 

n (l+i^i) = / exp(-e®) e® ds (3.25>’ 

— OO 

This is v/hat exactly is the second integral in eqn, (3*22). 
Fortunately values of Gamma function are available in 
tabulated form. Defining K(p) as 

K(p) - P(l+ip.) 

one gets from eqn. (3,22), 

F(p) - K(p,) . G(p) 
or 

G(vO = (3.26) 

If one takes the inverse Fourier transform of eqn, (3,26), 
then one gets 
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g(«-’') = 7 G(|i) * (3.27) 

•-OO 

Hence a plot of g{e~^) with respect to y will give the value 
of exponents as peaks, the amplitudes of which are proportional 
to N^/a^* 

3.3,3 Implementation and Results: 

The above algorithm is implemented using BASIC 
program on the signal analyser (IWATSU model SM2100B). The 
integrations are done by Simpson’s method. The listing of 
the program is given in Appendix II, To prove the method, 
both analytical and experimental data were taken. Figs, 3,2 
to 3,5 illustrate the results of the analysis,. The analysis 
is in the form of a plot of g(a} as a function of a. However 
such a plot cannot be directly obtained in the signal 
analyser, primarily because of the fact that a given function 
can be displayed only at discrete instants of time. The 
functional values at each of these instants are stored in 
the form of an array in the m^aory. These array locaticais are 
accessible through software. The above fact is exploited by 
storing a and g(a) values simultaneously in two different arrays^ 
And when the two arrays are displayed, the values read by the 
cursor at each instant give the values of, g(a) and a. 

To appreciate the multiexponential analysis, four 
examples in Figs. 3,2 to 3.5 are used to illustrate its 
various details; 
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1* All the plots jjn. each figure show a spectra like 
behaviour* This is so because of finite cut offs for x 
and p variables in e<jn, (3,18) and eqn, (3,27). For instance 
eqn, (3,18) with finite cut offs for x and p, is given by; 

t^o 

F(lJ>) » / . e^f(e^) exp(ipx)dx + E[x^,p] (3,28) 

-^0 

The extra term above is an indicator of the error due to 
cut offs in X and p. This error manifests itself in eqn. 

(3,26) where F(p) is divided by K(p), As the values of 
p, are increased, the K(p,) values decrease very rapidly and a 
small error E(x^,p) in F(p,) gets anplified to large error 
peak. This can be seen in Figs, 3,2 to 3,4 where as p, is 
increased from p=3 to p-5 and then to p=7, the number of 
peaks increases. In Eig, 3.5, the same thing is observed 
for p=2ii p=4 and p~6, 

2, ’The existence of error peaks in these can lead to 
loss of peaks, corresponding to actual values. However 
the redeeming feature of this method is that as one varies 
values of p, the error peaks shift their position while real 
peaks do not. This is shown in each of the figures and for 
each of the p* s by placing the cursor at the peak or peaks, 
which move very little. 
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3, The resolution of each actual peak is neximum at 

higher p*s and hence accurate ® known by increasing 

p* s gradually so that these are properly distinguished at 
lower p*s when error peaks are smallest and accurately 

known at higher p's, where their resolution is maximum, 

4, The variation in a in all figures 3,2 to 3,5 is 

from 0 to 1, Hence all decay constants need to be scaled 
in time so that a fails between 0 and 1. For instance in 
Fig, 3,4 an exponential of the form lOOe*"^^^^^ is used. For 
decay constant to fall between 0 and 1, it is scaled in 
time as The result of the analysis then 

would give the peak of g(a) at a« ,04, The value as read 
from Fig, 3,4(c) shows this a value as .040649, Similarly 
scaling is done in Fig, 3.5, In both the figures hence the 
scaling factor used is lOOOOO, As the exponentials of 

Fig, 3,2 and 3,3 have exponents less than one, therefore 
no scaling is needed in these cases, 

5, For the case wSiere more than one exponent is present, 
should increase as increases. For our illustration, 
the two double exponentials arrf ^ 

2e~^00000t 

An Another precaution which needs to be taken is that 
the given exponential function should be properly sampled. If 
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the whole exponential infoxn^tion is not available as input 
to the algorithm, error is found in the determined values. 
Fig, 3,6 shows how this particular effect in case of two 
different cut offs on the exponential gives a different a 
value. For the transient lOOe”"^*^^^, the values of e^f(e^) 
are plotted in (a) and (b).. The values along the abscissa 
are proportional to x values. In 3, 6(a) the sampling is 
done till the exponential has totally decayed and in 3, 6(b), the 
data obtained is only upto an earlier time instant, where the 
transient has not decayed completely. The results of the 
analysis show that for the case in 3.6(a), the peak occurs at 
4064,9, This is seen in Fig, 3, 6(c) while Fig, 3,6(d) shc^ws 
a value as 5493*1 for abrupt cut off case, 

7* Finally an the values found by multiexponential 
analysis in Figs, 3,2 to 3,5 are compared with the actual 
values. In Fig. 3.2, for analytical function lOOe ' 
the analysis is done at p=s3, p=5 and p=7, CMly one peak is 
seen to move very little alcaag the abscissa. This peak 
is identified by the cursor and its value is read in (c) 
as ,020141. 

Fig, 3,3 analyses double analytic exponential 
le~*^^^^ + The p value is varied from p=3 to 

p=5 to p=7. There are two peaks found to move very little 
and they are read in (c) and (d) as ,018 and ,082 respectively. 
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In Fig, 3,4 experimentally obtained single exponential 
is analysed at p=3, p.=5 and pis7. One peak as identified by 
cursor is shown to be almost stationary. The value of a at 
this peak is read in Oc) as ,040642. As the scaling factor 
of 100000 was used, the actual a value is 4064,2. 

Lastly the double exponent and 

in Fig, 3,5 is chosen particularly since it v^as found using 
conventional methods that this was numerically equivalent 
to its least square fitted expression 2.98375 e*"^056,71t 
shown in Table 3..3, Now as the multiexponential analysis 
is done on it, the two exponents show themselves clearly 
as p, is varied fron 2 to 6 as shown in Fig, 3, 5(a) and (b) 
respectively. At p=2 there is seen a broad peak, which 
breaks into two as p is increased to 4 and then to 6, The 
values of a^’s at these peak points are read jn Fig, 3, 5(c) 
and (d) as ,054931 and 1,6528. Using the fact that 
scaling factor of 100000 was used, the actual values 
are 5493 and 16528. 

Thus it is seen that multi exponential analysis is 
very convenient to detect inui tie xponents since the number 
and values of components fall out of the analysis. The 
fallowing subsections discusses its application to 
relaxation data. 
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3.3.4 Application to Relaxation Data: 

Application to particular case of isothermal relaxation 
is considered here. The following points can be made, 

1. This method is used to ascertain the number of 
components in the transient and get their approximate values, 

2* Because of the limitation listed in (6) in Section 
3.3.3, the results from an unknown sample need to be 
considered carefully, 

3. p, values need to be varied from the lowest possible 
to the largest possible values so that the true peaks are 
unambiguously distinguished* 

4, Major problems faced when implementing the algorithm 
on the signal analyser are its long execution time and 
limited accuracy. Both Imitations ccaitribute to the 
determination of only approximate values, taken over rela- 
tively long times. Needless to say given a faster and more 
accurate machine, it is possible to get fairly accurate 
spectra, 

3.4 ANALYSIS FOR K.TS; 

IXTS, as explained in Oiapter 2, is an elegant 
procedure for characterizing deep levels, even in the 
presence of multiple traps, which have reasonably apart 
activation energies. It has advantages over the above 
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analysis procedures, since it is spectroscopic in nature. 

At the same time it retains the essential features of 
space-charge relaxation methods. The conventional DLTS, as 
mentioned in Chapter 2, uses inefficiently the data and 
requires a double Box-car integrator. For each point in 
Arrhenius plot, a separate thermal scan is needed. The time 
and effort needed for this are prohibitive for very large 
tenperature variations. This disadvantage is overcome by 
separating the acquisition** of signal from its analysis [8], 
In this modified procedure, the isothermal capacitance and 
current are acquired and stored in memory of signal analyser. 
Corresponding to setting of rate window in Box-car integrator, 
a program is written, wherein variation of a single variable 
changes the rate window, A ELTS spectrum is thus obtained. 
However me should be careful in interpreting this data, 
since now temperature is held constant and not varied 
continuously. 



CHAPTER 4 


M£ASUREi/£NT SYSTEM FOR DA.TA ACqUISITION 
4,1 INTRODUCTION; 

This chapter deals with the measurement aspects of 
capacitance and current transients. Both the signals have 
very low amplitude and hence are very noise prone. The 
capacitance measurement, particularly, has to be done at 
500 KHz or above, since in order to monitor space charge 
capacitance, the monitoring signal should not disturb the 
transition capacitance. Section 4,2 deals with the measu- 
rement system used for data acquisition. Section 4,3 and 
4,4 discuss, in particular, the design of capacitance and 
current measurement modules-, 

4., 2 MEASUREMENT SYSTEM FOR DATA ACCIJISITION ; 

The measurement system can be divided into four 
blocks, as shown in Fig, 4,1, A brief discussion rs given 
for each block, 

1, Cryostat; 

This essentially contains a Liquid Nitrogen dewar ■ 
and a sample holder to place the sample at various depths 
in the cryocan The tenperatuxe variaticari can be done 
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exploiting the temperature gradient within the devvar. The 
temperatures above room temperatures can be obtained with 
the help of a soldering iron, 

2, Temperature Measurement Block: 

Temperature measurement is done with the help of a 
fine quality copper.- Cons tan tan thermocouple. The 
thermocouple junction temperature is monitored directly 
with 1iie help of a Digital multimeter (Keithley model 191). 
Room temperature correction is done after a mandatory warm 
up time of an hour or so. This is done by nulling this 
initial reading on the multimeter itself. An error of 2 to 
3®C can result in absence of this correction, 

3, Signal Acquisition Block: 

The acquisition block consists of following acce- 
ssories, each of which is briefly described. 

D.C, bias supply : This is needed in case of 
minority carrier pulse when the sample needs to be slightly 
forward .biased. It is normally superimposed on the pulse 
bias. 

Pulse bias supply : Tnis is used to obtain the 
majority and minority carrier pulses. Here the pulse should 
be negative going from 0 V to -10 V, The pulse generator 
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used should have pulse rise and fall tiraes of the orders 
of microseconds. In this case, Yamuna Pulse Generator 
has been used* 

High frequency slnevvay e Generator : This is 
particularly used in case of capacitance transient measu- 
rements, The exact need for it is explained in the design 
of capacitance measurement module. The model used is 
GR-1210C. 

^sp^citance measuremen t module : This is used 
for getting both majority carrier and minority carrier 
trap transients. High frequency (500 KHz or more) sinewave 
is used to monitor the capacitance changes in the space 
charge region. Its design is explained later on in Sec,4,3, 

Current measurement module ; It is also used to get 
transients from majority and minority carrier traps. Its 
design details, too, are explained later on in Section 4,4, 

4, Recording and analysis block: 

These functions are performed by IVfATSU signal 
analyser (model SM-2100B), The options necessary for 
our applications are: 

Averaging ; It is needed to improve the S/N 
ratio of very noisy capacitance and current transients from 
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the sample diode. Summation averaging is used in this 
case, 

Basic option; This option is used to obtain the 
various analys s programs described in chapter three. 

Minifloppy option ; This serves as our secondary 
storage unit for the current and capacitance transients. 
This helps in separating the acquisition from analysis, a 
desirable characteristic for DLTS, 

4.3 DESIGN OF CAPACITANCE TRANSIENT MEASUREMENT iVDDULE: 

The aim of this module is to measure the capacitance 
change in the transition region of the sample, during the 
relaxation process. At the outset, one realizes that the 
method used to measure capacitance should not disturb the 
charge condition in the depletion region. Hence the 
monitoring signal should be both s..:all in magnitude as well 
as should be of high frequency such that charge change 
should be slower than it. Here this transient is measured 
using principle of phase-rsensitive detection, 

4.3,1 Theory of Phase Sensitive Detection: 

The principle of phase sensitive detection can be 
explained using Fig, 4,2. 

Here there are two channels one is the reference 
channel and other is the channel through system under test. 




Fig. 4,2; Block diagram of a phase sensitive detec-tor 






Both the channels are locked-in to the same carrier signal 
frequency Hence the name lock-in detection also 

applies to the method. Calling the input signals to the 
multiplier as and X^, then 


cos oj^t 

^2 = ^ cosCwct -h 0) 


Output of the multiplier win be X^ = X 2 = 


cos(u^t + 0 ). 


E E 

Xq = -^JS [cos 0 + cos i2b)^+ 0)] 


(4.1) 


If we filter the output of the multiplier such that << 

2w^ then 


>^0 = 


ErEm 
c m 


cos 0 


(4.2) 


If in eqn. (4.2) noWj a phase shifter is introduced in the 
reference channel such that phase difference between the 
two signals and X 2 is zero, then 


V 


E E 
c*'m 


(4.3) 


Hence if the ‘system under test’ is so designed that changes 
in capacitance are reflected as amplitude modulation of E^^^, 

X the output signal v/iil be proportional to capacitance 


then 
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changes directly and vi/e have succeeded in getting capa- 
citance transients, 

4,3,2 Design of Setup: 

Keeping the above block diagram of Fig. 4,2 in mind, 
the capacitance measurement module should look like as in 
Fig, 4,3* A brief discussion of each block is given along 
with the realization. 

Summer ; Summer is needed to sum up three 
quantities; 

1. Pulse bias 2. D,C. bias (for minority carrier 
pulse) 3, High frequency sinev^^ave. 

As operation needs to be done at high frequencies, 
LM318/CA3140 were found suitable for this application, 

I/V converter ; The design of this is the trick iest, 
since the. signal level is quite low here and this being the 
first stage following Fritz's law, it needs to be as noise 
free as possible. It is instructive to know the order of 
currents available to be amplified by this converter. 

Amplitude of carrier signal = lO mV 

Frequency of carrier signal = 500 KHz 

Typical capacitance of 1N4148 

diode (from data sheets) = 4 pF 



Surmner Sample diode j/Y converter Low noise 
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Capacitive reactance of the diode 

1 

271x500x1^x4x10“^^ 

= 80K 

-3 

Current through the diode = ^22^:0—^ 

80x10'^ 

= 100 nA 

This current has to be converted to voltage at 500 KHz and 
further amplified. The noise itself can be in millivolts 
and hence one needs to be careful. Here this is obtained 
by following steps: 

a) This low current is dropped across a precision 
resistor. The value of this resistor is so chosen that it 
is much Smaller than the impedance of diode in reverse bias 
and at that frequency, 

b) The developed voltage is then buffered using a 
very low noise and excellent high frequency voltage follower 
LM310. 

c) The buffered voltage is then amplified using a 
wideband amplifier (HP model 46lA), having excellent input 
characteristics. However only gain of lO is possible 
because amplifier oscillates due to mismatches at output. 
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d) Tns signal level still is low and hence Idie 
oscilloscope vertical amplifier output itself is taken. 

This output is applied to the multiplier block as one 
of its inpu ts . 

Phase shifter : It is achieved using CA3140 in a 
phase-shifter configuration and is detailed in the actual 
circuit diagram, 

D) Multiplier ; The process of multiplication can be 
done using either a four quadrant multiplier or a balanced 
modulator. However the disadvantage with the available 
balanced modulators is that output signal level is small 
and hence further amplification is needed. Hence a fou3>- 
quadrant multiplier XR--2208 is used in its phase- sensitive 
detection mode as shown in Fig, 4.4, 

V0 in this case is given by; 

V0 » cos 0 

where is the phase detector conversion gain. If 
> 50 mV r.m.,,s., = 2V/radian. However, if input 

signals are <50 mV rms., then changes as the input changes. 
However, the ouiput from Fig. 4.4 is differential in nature 
and hence a difference amplifier must follow. Again LM318 
Was found quite satisfactory for this purpose. 
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Since is small, its amplitude is increased by 
amplifying the reference signal in reference channel since 
any way on that channel, signal has to be larger than 50 mV 
r.m.s. The signal circuit then looks like as in Fig. 4.5, 

4.4 DES3BN OF CURRENT TR/sNSIENT MEASUREMENT MODULE: 

Current transient measurement module has tv/o features 
different from that of capacitance measurement module in 
previous section. One is that one need not work at high 
frequencies in iiiis case and second that level of signals 
is smaller in this case than of capacitance transients. 

The latter is particularly true since now it is the reverse 
bias current, which has to be monitored. In silicon diodes, 
this current is of the order of a few nanoamperes. Again 
due to requirements of Fritz's law, the first current to 
voltage converter should be as noise free as possible. The 
block diagram of such a setup is in Fig, 4,6, A brief 
description of each block is provided, 

a) I/V Converter ; This stage should be as noise free as 
possible. Hence a FET input low noise OPAMP is used, 

CA3140 is characterized by low offset voltage, low offset 
current high CMSRR and high input impedance. Several stages 
of tran«impedance are provided. 
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VO-X-'to qe 3 lapl if ie r ; Again an amplifier with low noise, 
low offset current and voltages should be used. Multiple 
gains are provided in this case, too* 

The final current transient circuit is as shown in 
Fig. 4,7, Some very important precautions need to be taken 
regarding: 

a) layout and connections 

b) selectim of passive components 

For layout, the above low. noise, low level current transient 
set-up or electrometer should have proper shielding and 
guarc'rings must be provided to avoid loss of signal [13,14], 
Both capacitance and current transient signal measurement 
Setups have been developed and used to acquire these trans- 
ients from a p"^ junction diode. The details of this ^re 
in Chapter 5,- 



CHAPTER 5 


EXPERIMENTAL RESULTS AND DISCUSSICN 

5.1 INTRODUCTION .* 

In this chapter the details of current and capacitance 
transients from commercially available Gold doped P'^N junction 
diode have been given. These transients have been analysed 
to obtain the values of activation energy and capture cross- 
section related to Gold acceptor in silicon. The details 
are described below: 

5.2 CURRENT TRANSIENTS: 

The experimental system used to measure current 

transients has already been described in Fig, 4,1, A total 

of eighteen such transients spanning a temperature range 

of 110^ from 236\'to 274^^ are taken and those.. at two temporal 
are shown in Fig. 5,1 and 

/waveshapes are similar to the theoretically predicted one 
in Fig. 2*4. Before analysing these for determination of 
activation energy and capture cross-section, following 
considerations, over and above those already described for the 
measurement system of Fig, 4,1 are taken. 

1, Each of the transients is measured using a majority 
carrier pulse wherein the voltage switches from O V to a reverse 
voltage -Vo, To illustrate the shape of pulse bias, the 
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majority carrier pulses in case of transients at T = 251^ 
and T =* 274®ic are shown in Fig. 5.1, 

2. It has been found that pulse bias at a single 
frequency cannot be used for the complete temperature range 
used* This is because at lower tenperatures , the time of 
decay of current transients are so large that it does not 
totally become zero, leading thus to conditions for DDLTS 
(Double DLTS). In such a case the analysis needs to be 
changed accordingly, Qi the other hand at high temperatures, 
the decay is so fast that a set of indistinguishable 
transients are obtained close to the negative going transi- 
tion of bias pulse. Therefore for higher temperatures, high 
repetitive rate pulses are used and for lower temperatures, 
the repetitive rate is made so large that the decay is 
complete. This particular point is illustrated in Fig, 5,1 
where the transient at T=274\ is measured using a puise of 
about 2,40 ms period and one at 1=251*^ obtained with 
about 85 ms pulse period. 

The transients obtained are submerged in noise. Hence 
the signal averaging capability of the signal analyser is 
u»ed to increase SNR for these transients and it has been 
found that at least five hundred such averagings are 
necessary. 
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Before analysing the data, it is necessary to check 
out some features of current transients which are predicted 
by the eqn. (2.26), These are: 

a) that the decay constants of the current transients 
decrease with decreasing temperature and 

b) that the initial anplitude becomes smaller as the 
temperature is lowered. 

Property (a) above is verified in Fig. 5.2 where the 
current transient signals at three different temperatures 
(T=346°K, 325*^ and T=306°K) have been plotted. In Fig, 

5, 2(a), T=325*\ and 1=346^ transients are compared and 
values are read by cursor at t = 261.23 ps to show a 
higher value of 826.43 mV for T=346°K transient meaning 
thereby that lower amplitude for latter transient is due 
to its faster decay. Similarly Figs. 5,2(b), 5.3(a) and (b) 
illustrate this property. 

Theproperty (b) above can be observed by comparing 
the order of magnitudes of transiaits in Fig. 5,2 with 
those in Fig. 5.3, For instance at the cursor setting of 
t=26l,23 p,s in Fig, 5,2, the ma^itudes are in hundreds of 
millivolts while in case of Fig. 5.3 at t=S.oa78 ms setting 
of cursor, the magnitudes are in terms of millivolts. The 
temperatures in Fig. 5.3 are lower than temperatures in 
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Fig# 5#2, meaning hence that amplitudes decrease with 
temperature, 

5.2 ANALYSIS FOR CURRENT TRANSIENTS: 

The eighteen transient signals experi- 

mentally obtained, have been analysed for determination 
of activation energy and capture cross-section. Many of 
•the methods [4], [5], [8], [10] including the DLTS assume 
that the transients are single exponential. In practice 
it has been found that sometimes multiexponentials occur. 
Hence efforts should be made to find whether the transients 
are single exponential or not. The me1±iod described in 
Chapter 3 for multiexponential analysis is used for the 
first time to get this preliminary infornetion. In case 
the signal shows single exponential behaviour, then either 
DLTS or any of the other single exponential analyses can be 
used. If it is multiexponential, then this algorithm itself 
can be used to get Arrhenius plot. 

Therefore our analysis of current transients starts 
with multiexponential analysis of these transients. The 
details of this are given below: 

a. Multiexponential analysis: 

Herein the analysis is done at two temperatures. The. 
current transients at 321.63°K and 258.5°K are chosen for 
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specifically this effect is illustrated in Fig. 5.6(a) 
and (b) at tw© different temperatures. 

2. Because of the above raeitioned problem, using 
arguments of chapter three, area method is used to find 
values of B, These values for each temperature are given 
in Table 5,1, As can be seen, the values of B decrease 
with lowering of temperature, because it is due to 
leakage current which is temperature dependent, 

3, The values of a and A are obtained from iiie original 

curve minus the value of B (i.e, Ae” The program for 

this is provided in Appendix I, The values of A and a are 
also tabulated in Table 5,1, The values of both decrease 

with decreasing temperature according to egn, (2.26), 

* 

4. These values of a with corresponding temperatures 

T are used to obtain the Arrhenius plot. As motioned in 
Chapter 2, the slope of this plot gives the activation 
energy and intercept the capture cross-section. The values 
as obtained from the graph are shown therein, 

5, The above values of a and l/T are also fitted by a 
linear regression algorithm provided in Appendix I, The 
values of E, the activation energy and g , the capture 
cros s— sec tion are shown in Table 5',2# 

Next we discuss the capacitance cfctat 
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Table 5,1: Current Transient Data 
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5.3 CAPACITANCE TRANSIHnITSs 

There is only one change while measuring the capa- 
citance transients as compared to that of current transients. 

In the measurement system of Fig, 4,1, the current measurement 
set up is replaced by capacitance set up, A total of eleven 
measuroments axe made over a 50°K range from 298,9^ to, 

250*^. The waveshape of typical capacitance transients is 
illustrated in Figs, 5,7 and 5,8, These waveshapes are found 
to be same as theoretically expected in Fig, 2,4, Both the 
considerations mentioned fox current transients hold in this 
case too. This means that we have again majority carrier 
pulses, switching from 0 V to a reverse voltage - Vj^. This 
is illustrated in Figs, 5,7 and 5,8, vhere the pulse bias 
is shown below the capacitive transients at four different 
temperatures. Due to the same difficulties as faced in 
current transients, a set of different repetition rates are use 
For tho highest temperature of those at T=292,71^ the pulse 
period is about 9.0 ms and for the lowest temperature 
T = 246^ illustrated in Fig. 5.8(b), the pulse period is 
90 ms. The intermediate two temperatures similarly have : 

pulse periods of 16 ms and 45 ms, 

. . ■ f 

The capacitance transients are again very noisy and 
at each of the temperatures imilt ip oint averaging is done to 
get cleaner signal. The effect of noise is more at lower 
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temperatures when signal strength becomes small. This 
effect can be observed in Figs, 5.7(b)^ 5.8(a) and (b) 
and Fig, 5, 9(b), 

In this case too, we check the capacitance transients 
for some features predicted by eqn, (2,19), These are: 

a. that the decay time increases as the temperature 
is decreased 

and b, that the initial amplitude remains constant. 

Property (a) is verified in Fig, 5, 8(a) and (b) and 
Fig. 5, 9(a), In Fig, 5,8(a), the capacitive transients at 
T = 298,8^ and T =: 288*43®K are compared with each other. 

At the setting of the cursor at t = 830,07 ps, the values 
read are 88,195 mV and 91,735 mV respectively, meaning 
thereby that the transient at T=298,87°K decays faster than 
at T=288,43°K, The same effect can be seen at other 
comparisons shown in Figs, 5, 8(b) and 5, 9(a), 

property (b) can again be verified in Fig, 5, 8(a) 
and (b), where initial amplitudes are found equal, 

Vte proceed next with the analysis fox deter- 
mination of activation energy and capture cross— sections, 

5.3.2 Analysis for Capacitive Transients: 

Here again as in current transients, a serious need is 
fgH; -Iq verify if the capacitive transients are single- 
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exponential or not so tl^t either ELTS or any of the other 
single exponential analysis could be used or not. Therefore 
the multicxponential analysis is done to verify the existence 
of the single exponencity. In case the transients are 
single exponential, then either DLTS or any of the other 
single exponential analyses can be used. If, however, 
these are multiexponential, then this algorithm can be itself 
used to get Arrhenius plot. 

Therefore our analysis of capacitance transients starts 
with multlexponential analysis. The details of this are 
given below: 

a. Multicxponential analysis: 


■ -The, capacitance transients at two temperatures are 
analysed for multiexponencity. These transient are at 
T=252.53^K and T=:2'r9.6l'^ and are shown in Figs. 5.10(d) 
and 5.11(c) respectively. The multiexponencial Ig 
when applied to transient at T = 279. 61 °K gives for 
p=2, p=4 and p.6 in (a) and (b) and its value can be.e^d 
in Fig. 5.10(0) as 148^ This the transient at T=279.61°K 

has single exponential* 

Similarly the capacitance transient at T=252.53 “k 


is analysed for multiexponents, 

shown in Fig. 5.11 vAiere (a) and 


The results of this are 
(b) give thog(a) values for 



me 



Fig, 5,10: Multiexponential analysis cn a capac 





90 


p,-2» iJi-4 and [i=:6. One particular peak is seen to remain 
constant and the a value at this is read at the cursor point 
in Fig, 5,ll(c), The value found for a is 135.49, Hence 
the capacitive transient at 1=252.53^ is also single 
exponential. 

The reason why only two temperatures are analysed for 
multicomponents is tiiat the multiexponential analysis in its 
present implementation takes long time. Having verified that 
the transients are single exponential, any of the simpler 
and faster single exponential analyses can be used to find 
the values of decay constants. It is described below, 

b. Single Experiential analysis: 

Here the capacitive transients are assumed to have the 
form Ae‘““’^ + B, The following steps are taken to get values 
of A, a and B, 

1, The Area method is again used to find values of B, 
These values of Bare tabulated for a few temperatures in 
Table 5,3, As can be seen that this value is essentially 
constant because the capacitance at steady state reverse 
bias is same at all temperatures. These values of B are 
removed frem the whole waveforms at each of the eleven 
temperatures, 

2, The modified waveforms now obtained are of the form 
Ae“®’‘^ and values of a and A are determined using ratio of the 



Table 5,3; Capacitance Transient Data 
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exponential function at two instants. These values of A 
and a are also given in Table 5*3, The value of A does 
shew small change and it is both because of quantization 
erroi as well as the actual high frequency cut off in the 
signal analyser. The a values show a decrease, with 
decrease in temperature, 

3, These values of a with corresponding temperatures 
T are used to obtain the Arrhoiius plot. The slope of -ttiis 
plot gives the activation energy of the accept car level of 
Gold in silicon and is given in the plot itself. Though the 
intercept of this plot along the ordinate gives capture 
cross-section, this has not been got on the Arrhenius plot 
itself since the plot, then extends outside the graph, 

4, Instead values of the capture cross-section and 
the activation energy are accurately computed by linear 
regression algorithm. These results are tabulated in 
Table 5,2* 

A Comparison of these measured values with the ones 
published in literature is given in the next chapter. 



CHAPTER 6 
CGNaUSICN 

The values of the activation energy AE and capture 
cross-section g of gold acceptor level in gold doped 
silicon P'^N Junction diode 1N4148 have been measured using 
majority carrier pulses. The values of AE using isothermal 
current transient study (ITS) and using isothermal capacitance 
transient study (CTS) have been found to be '^0,53 eV and 
0.58 eV respectively. These values are well within the 
published values which vary from 0.52 to 0.58 [lO], Similarly 
the values of g found by ITS and CTS are ■X 0.8845x10*”^^ cm^ 
and 0,664,321x10”"^^ respectively , These g values are 
found well within those published and which vary from 9.0xl0”^'^cm^ 
to 2.5x10"”^^ cm^ [16] 

Current and capacitance transients Vi/ore obtained 
using especially designed measurement set ups capable of 
acquiring low level current and capacitance signals. 

A multiexponential analysis algorithm was used to help 
detect the multiexponents in relaxation data. 

Some limitations of the above work and ways to improve 
these are as follows: 

The multiexponential analysis, as implemented in our 
case, realizes Fourier transforms through numerical integra- 
tions, This means that an increase in execution time as well 
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as in numerical errors, if the FH routine of the signal 
analyser can be used, then both the above problems can be 
solved. 

An error analysis is particularly desirable for the 

multiexponential analysis. This is because the signal from 

the form Z A. e"*®i^ + B is reduced to S usinq 

2^1 i=l ^ 

Area method. The effect of such a step needs to be determined. 

The multiexponential analysis has the limitation that 
Aj^ should decrease with a decrease in a^. Due to laig 
execution time of the program, the extent to which this 
condition is stringent could not be known. V^ith FFT imple- 
mentation, it would be easier then to determine the extent 
of this- limitation. 

The current transient setup could be redesigned to 
remove the large current, component. The design is, however. 
Slightly difficult since this fast transient decays is loss 
than a few microseconds. 

Finally the capacitance module has been found noisy, 

A better layout and better choice of less noisy canponents 
can help solve this problem. 



APPENDIX I 

I, A PROGRAM TO FIND EMISSION RATE CCNSTANTS BY LINEAR 
REGRESSION ALGORITHM 

5 DIM X(.500), Y(500) 

7 H =: 2.44l4Er.06 
9 E3=50, E4=50 

10 J= = 55, J7 == 55 
12 T8 = 300 
15 T7 =3 30 
17 FOR J * 1 TO T8 
20 X(J) = H*(J1-E3) 

30 T9 =* B(2,0,J1) 

40 Y(J) a LN(T9) 

50 J1 = Jl+1 
60 NEXT J 

65 SI 0, Tl=6, S0«0, S2=0,T^0 
70 FOR 1=1 TO T7 
80 SI = SI + X(I) 

90 ■ 'tl = T1 +-¥<I) 

95 SO = SO + X(I)* Y{I) 

100 S2 = S2+X(I)*X(I) 

no T2 = T2+Y(I)*Y(I) 

120 NEXT I 

130 XI = S1/T7, Yl=Tl/n, X2=S2/T7 

135 Y2?=T2/T7, P = S0/T7 

140 P1=X2-(X1*X1), P2;= Y2-«(Y1*Y ) 

150 M = (P-X1*Y1)/P1 

160 C=Y1-M*XI 

165 PRINT *£3' , E3 

170 PRINT ‘MV, M, ‘O', EXP (C) 

172 T4 = 0 

175 FOR L=1 TO T4 
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180 Y(L) » M*X(L)+C 
190 IX =« EXP(Y(L)) 

200 12 = B(2,0,J7) 

205 G = (12-11) 

210 T4 + T4fG*G 

215 J7 « J7+1 

220 NEXT L 

230 R4 » T4/T8 

240 PRINT »MSE‘ , R4 

250 E3 « E3t-1 

260 IF E3 < = E4 GO TO 10 

270 STCS> 

T B PBOGBAM TO FIND EMISSION RATE OCNSTANTS BY FINDING RATIO 
AT TWO TIME INSTANTS 

5 DIM D(20) 

10 El a ? 

20 TO a 2.4414E - 06 
30 K a 30 
40 XkEI 
42. I « 1 
■ 45 Y a X+K 

60 TL = X* TO, T2, Y*TO 
70 Y1 = B(2,0,X), Y2pB(2,0,Y) 

80 XO = LN(Yl/Y2)/(T2-Tl) 

82 D(I) = XO 
85 X a Xi-1 
90 PRINT I, XO 
100 I a I+l 
no IF I < -4 goto 45 

185 S a 0 

190 FOR Jal TO I 

210 S a D( J) + S 

215 P a S/6 
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220 PRINT P,S 
230 STCa? 

I C LEAST SqUARE FIT PROGRAM TO RELAXATION DATA FOR 
ARRHENIUS PLOT 


5 DIM X(20), Y(20) 

10 T1 = 50.427, T2 = 44.95, T3 = 39.859, T4=33.6026, 
T5 =c 27.8684 

20 T6 » 20.3684, T7 = 14.04, T8=8.04, 19=4,06, 

Q1 =: -2.16, Q2 = -6.90 
30 Q3 = -12.186 

35 El = 9881.65, Efe7125, E3=5378 , 260, £4=3653.9, 

E5 == 2397 • 92 

40 £6 = 1400.6834, E7=704.34, E8=384.03, E9=269.417, 

Z1 » 179.109 

43 Z2 = 125.062,23=84.9573 
45 K=1.0E-6 

Z X(2)= 1/CT2+T), XC3)=1/(T3H) 

X(4) = 1/(T4+T) 

55 X(5) = l/(T5fT), X(6)=1/(T6)-I) 

57 X(7) = l/Cn+T),X(8)=l/(T^) i/roa-T), 

59 X(9) = 1/(T9+T), X(10)=l/(Q1+T)> x(ll) /(Q 

X(12)= l/(Cl3fT) vniWN(Z2') 

60 Y(9) = LN(E9), Y(lO)=LN(Zl), Y( 

Y(12)= LN(Z3) 

61 Y(l) = LN(El), Y(2)=LN(E2) 

62 Y(3) = LN(E3), Y(4)=LN(E4) 

63 Y(5) = LN(E5), Y(6)=LN(E6) 

64 Y(7) = LN(E7), Y(8)=LN(E8) 

65 S1=0, T1=0,SOO, S2=0,T2=0 

67 Ssl2 

70 FCR I==I TO 12 
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80 Sl=sSl + X{I) 

90 11=11 + Y(I) 

95 SO = SOfX(l)*Y(l) 

100, S2 = S2fX(l)*X(I) 

110 T2 = T2fT2+Y(l)*Y(l) 

120 NEXT I 

130 XI = Sl/S, Y1 =Tl/S,.X2=S2/S 

135 Y2 = T2/S#. P=s SO/S 

140 Pi = X2-(X1*X1), P2=Y2-(Yim) 

150 M = (P-X1*Y1)/P1 

160 C = Y1-M*X1 

170 PRINT 'M', M,'C» ,C 

175 D=M*8,61E^5 

177 V = (EXP(C))/(2.6E+26) 

180 PRINT ‘D‘ , D, ’V’ , V 
190 STC«> 

I.D PROGRAM TO FIND EMISSION RATE CONSTANTS USING RATIO 
METHOD 

3 El = 56, E2=74, E3=98 
5 TO = 2.4414B-6 
• 7 Tl-iTO^El, T2=T0*E2, T3sT0*E3 

10 K«8 V 

12 Y1=B{K,0,E1), Ya-B(K,0,E2). Y3=B(K,0,E3) 

20 XO = 1000.0 
30 PI = (EXP(-XO*Tl)) 

40 P2 = (EXP(-X0*T2)) 

50 P3 « (EXP(^X0*T2)) 

70 F1*T2*P2*'(P1-P3) 

80 F2 SB Tl*Pl*(P3-P2) 

90 F3 ss T3-*P3*(P2-P1) 

100 D1 « CPS-P2)*(P3-P2)v 
110 D2 = 



100 


120 E2 = ((P2-P1)/(P3-P2)) 

130 m = ((Y2-yi)/(Y3-Y2) 

140 N2 = E2-N1 
145 R2 = N2/D2 
150 XI ;= (XDU.R2) 

152 PRINT XI, N2 
155 V=ABS(X0-X1) 

160 IF V < l.OB-3 GO TO 190 

170 XO =s XI 

180 GO TO 30 

190 R = XI 

200 PRINT R 

210 STOP 

I.E PROGRAM TO FIND EMISSICN RATE CONSTANTS BY AREA METHOD 


5 

El = 

: 583, E2 s 1003 

10 

Ms3 


12 

H s 

2.4414E-6 

15 

J sr 

El, P=E2-1, K=l, E=0, 

17 

FI s 

= B{M,0,a)+B(M,0,E2) 

25 

IF K 

:=(INT(K/2))*2 GO TO 50 

30 

a 

+ B(M,0,J) 

40 

GO TO 55 

50 

E = 

Et.B(M,0,J) 

55 

J S 

J+1, KsK+1 

60 

IF J 

< = P GO TO 25 


65 PRINT • E’ , E, » ‘ , 

67 PRINT VFl*, FI 
70 S = ( Fl+2*Ef4^ ) * (H/3,0) 

80 PRINT S 
90 K « ABS ( S) 

100 Mis B(M,0,E1) 
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105 M2 = B(M,0,E2) 

110 R3 a K/Ml 
120 XO a 0 
125 T1 a E2<^H-E1*H 
130 FI a EXP (-XOm) 

140 M3 a M1-M2 

145 P2 = M3*(1-X0*T1-P1) 

150 P3 a (1-P1)*(K-M1*T1)*X0 
155 R2 = {M3/(1-P1) 

162 V2 a Ml-Vl. 

163 PRIKT XO, R2,V2 
165 XO a X04-500 

170 IF XO < a l.OE+5 GO TO 130 
180 STOP 

I.F PROGRAM TO FIND EMISSION RATE CONSTANTS BY RATE 
WINDOWING CC3NGEPT 

10 DIM D(60,20), W(60) 

20 J a 1, D1=2.4414E-6, L=1 
25 J a 1, X=5, N=1.5, Nl=2 
60 FOR lal TO 8 
65 PR. J 
70 FREAD (I,J) 

80 J a J+l 
90 . NEXT I . 

100, £1=132, E3all9 

110 E4a364 

120 M=1 

130 E1=E1+X 

140 E2aN*{El-E3)+E3 

150 T1 a D1*(E1-E3): T2=D1*( E2-E3) 

160 T0a((Tl-T2)/LN(Tl/T2)) 

170 W(L) a l/TO 

180 FOR K=i TO 8 

190 PI a B(K,G,E1)-B(K,0,E2) 



200 D(L,M)=P1/B(K,0,E3) 

210 M =t Mfl 

220 NEXT K • 

230 L L+1 

235 PR. ‘NIS N1 

240 IF E2 < * E4 €0 TO 270 

250 N1 = Nl-l 

260 GO TO 280 

270 IF N1 > =2 GO TO 120 

275 IF Nl=l GO TO 355 

280 IF N1»0, GO TO 370 

290 FOR Ml=l 

295 PRINT J 

300 FREAD (M1,J) 

310 J=lJ+l : 

320 NEXT Ml 
330 L=1 

340 El«132, £3*119 

350 E4«364, Mss9 

360 GO TO 130 
370 S a 1 

380 PR. ’WINDOW, W(S) 

390 FOR N=1 TO 12 

391 S5aO 

392 FOR 32=1 TO 150 

394 35 a 35+1 

396 NEXT 32 
400 PR. N,D(S,N) 

420 NEXT N 
430 3 a Sfl 

440 IF S < a L GO TO 380 
450 STOP 
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PROGRAM TO FIND MULTI EXPONENTS IN A DECAY TRANSIENT FROM 

RELA)(ATION STUDY 


10 REM *******^t** 

20 REM PROG. FOR MULT I EXPONENTS 

30 REM ********^* 

40 DIM X(IOO), Y(IOO) 

45 DIM Z{225),T(225> 

50 X(l) = 1.000: X(2)=0. 9902145: X(3)=0. 9619537 
X(4) = 0.9182789. 

60 X(5) = 0.8633796: X(6) = 0.801 6901: X(7) = 0.7371 572 
X(8) = 0.6728157 

70 X(9) = 0.61081 : X(10) = 0.5523: X(ll) = 0.4980: 
X(12) = 0.4480. 

80 X(l3) = 0.4022: X(l4) = 0.36036: X(15)= 0.32210 
X(16) = 0.28713 

90 X(17) = 0.25513: X(l8)= 0.2258; X(19) = 0.19801 
X(20) = 0.17440 

100 X(2l) = 0.1519: x( 22 )= 0 . 13138: X(23) = 0.11272 
X( 24) = 0.0958 

110 X(25) = 0*08055: X(26) = 0.06687: X(27)=0. 05470 
X(28) = 0.043925 

120 X(29) = 0,03449: X(30) = 0.02631: X(3l) = 0.0193( 
X(32)= 0.01335. 
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130 X(33) = 0.00840: X(34)*0. 00434: X(35)=0. 00109 : 

X(36)= -0.001446. 

140 X(37) « -0.00335: X(38)= -0.004713: X(39)= -0.005609: 
X(40) = -0.0061. ; 

150 X(41) -0.006305: X(42)=: -0.006239; X(43)= -0.005978 

160 X(44) « -0.005572: X(45)== -0,00506: X(46)= -0.0045017: 
X(47) -0.0039079. 

170 X(48) a -0.0033124; X(49)=! -0.0027360: XC50)= -0.0021946 
X(51)= -0.0017. 

180 X(52) =s -0,001258: X(53)= -0.00087; X(54)= -0.00549: 
X(55) = -0.0002819. 

190 X(56) = -0.0000687: X(57}= -0000946: X(58)= 0.0002131: 
X(59) = 0.0029. 

200 X(60) = .0003392: X(6l) == .0003588; X(62)= .0003569: 
X(63} = .0003390. 

210 X(64) = .0003096; X( 65) "= .0002729: X( 66)= .0002323: 
X(67) = .0001906 

220 X(68) = .0001501; X(69) = .000112: 

230 X(70) = ,0000785: X(7l) = .0000493; X(72)= .0000250. 

X(73) = .0000055. 

240 X(74) = -.0000093! X(75)= -.0000199! X(76)= -.0000256. 

250 X(77) = -. 0000303 ! X{78)= -.0000313! X(79)= -.^.jC03C4 

260 X(80) = -.0000280! X(81)= -.0000247 


270 Y(l) 0.0000: Y(2)= -.0.0368213: Y(3)= -.0.108484: 
Y(4) = -0.150849 

280 Y(5) ss — .lol45: Y(6)as -.loo64: Y(7)= -.20619: 

Y(8) a -.20288. 

290 Y(9) =. -.19177: Y(10) = -.175144; Y(ll)= -.1549509: 

Y(12) a -.132851. 

300 Y(13) a -.110161; Y(14)= -.087884: Y(l5)s= -.0667373: 

Y(16) = -.047204. 

310 Y(17) = -.0295912: Y(18)= -.0140620; Y(19)= -.0006752. 

320 Y(20) = -.010583: Y(2l) = .019804: Y(22)= 0.027107 

Y(23) a .0326383. 

330 Y(24) a .0365665; Y(25)a .029067: Y(26)= 0.0403226: 

y(27) a 0.0405085. 

340 Y(28) a .0397995: Y( 29)= .0383604: Y{ 30 )a 0,0363451: 

Y(31) = .0338961. 

350 Y(32) a .0311420: Y(33)= ,0281976: Y(34)a 0.251 629: 

Y(35) a .0221241. 

360 Y(36) a .019153: Y(37)a .0163071: Y{38)a .036326: 

YC39) = .0111631 

370 Y(40) a ,0089212; Y(4l)a ,0069205; y( 42)= .0051654: 

Y(43) a ,0036545. 

380 Y( 44) a .0023793: Y(45)a ,0013268: Y( 46)= .0004807 : 
Y(47) a -.0001781. 
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390 Y(48) = -.00067: 7(49)= -,0010167 1 Y( 50 )= -.0012393: 

Y(51) = -.001358 

400 Y(52) = -.0013942: Y(53)= -.0013645: Y(54)= -.0012859. 

410 Y{55) = -.0011729: 7(56)= -.0010382; Y(57)= -.0008923 

420 7(58) = -.0007437; Y(59)= -.0005993: Y(60)= -.0004642 

430 Y(6l) = -.0003418: Y(62)= -.0002343 : 7(63)= -.0001429 

440 7(64) = -.0000678: Y( 65)= -.0000083: Y(66)= ,0000366 

450 Y('67) = .0000683 : 7(68 )= .0000887: Y( 69)= .0000996: 

7(70)' = .0001026. 

460 7(71) = .0000997; Y(72)= .0000925: Y(73)= .0000 822: 

7(74) = .0000703 

470 7(75) = .0000577: Y(76)= .0000452 

480 7(77) = ,0000335 : 7(78) = .0000230: Y( 79)= ■.0000139: 

7(80) = .0000063: 

490 7(81) = .0000004 

500 SI = ^.01: S2 = pa<S\ 

502 E2=5i 5: J2=1: X7=0: El=2.4414E-6 

503 E3eHt)5: E4=i65: E5=237: E6=393: E7=489: E8=681; 

175 . StiS 2m 

E9=1065 

X ^f\ 

509 FOR ^ = 1 TO ^ 

510 X7 = El * I 

511 X8 = lOOOOO* X7 

512 T(I) = LN(X8) 



513 PR. T(I), B(8,0,I+E2), J2,I 

514 Z(I) =r X8 * B(8,0,I+E2) 

515 J2 = 32+1 

517 NEXTI 

518 PR.I,J2 

520 X7 = El * J2 
525 X8 = lOOOOO* X7 
530 T(I) = LN(X8) 

532 PR.T(I), B(8,0,E2+J2),J2,I 
535 Z(I) = X8 * B(8,0,E2 + J2) 
540 1=1+1 

542 IF J2+E2 < E3 00 TO 580 
545 IF J2 + E2<E4 C30 TO 590 

547 IF J2+E2 < E5 GO TO 600 

548 IF J2+E2 < E6 GO TO 610 

550 IF J2+E2 < E7 GO TO 620 

560 IF J2+E2 < E8 GO TO 630 

565 IF J2+E2 < E9 GO TO 640 

570 GO TO 650 

580 J2 = J2+2: GO TO 520 
590 J2= J2+4 ; GO TO 520 

600 J2 = J2+6: GO TO 520 
610 J2 = J2+8; GO TO 520 
620 J2 = J2+16: GO TO 520 
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630 J2 = J2+32 : GO TO 520 

640 J2 = J2+64 ; GO TO 520 

650 PI = I-l 

660 PR. ‘Pl’» PI 

670 U = 0: Cl = 1: IA=4 

680 CHAIN 31 


File 31 ; 


10 

20 

30 

35 

37 

2 ? ^ 

40 

50 

60 

70 

80 

90 

91 

92 


REM CONTINUATICN FOR MIN PROGRAM 
REM 

dim d(225), E(225) 

DIMV(lOO), F(IOO), G(IOO), C(lOO) 

FOFa'''= 1 TO PI 

T3 = Cos(U*T(J)) : T4 = Sin<Ml'*T( J)) 
D(J) = T3*Z(J): E(J) = T4*Z(J) 

NE)Cr J 

B(2,0,C1) = U 
GOSUB 370 
S3 = S5 

FOR J9 = 1 TO PI 


94 D( J9) - E(J9) 
96 NEXT J9 
100 GOBUS 370 

105 34 = S5 
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110 C(G1) = S3-*^(l/(SQR(2*PI))) 

120 F(Cl) = S4*(1/(SQR (2*PI))) 

130 U = U + 0.1 
140 Cl = Cl+1 
150 IF U < SI GO TO 40 
160 FOR K=1 TO 7rS2- ( ^ 

170 D3 = X(K) : D4 = Y(K} 

180 12= ((X(K)) (X{K)}+(Y(K))*(Y(K))) 

190 G(K) = ((C(K}™) + (F(K}--r)4))/l2 

200 V(K) = ((F(K} D3) - (C(K)-‘W) )/l2 

210 NEXT K 

220 Y1 = 0 ; Y2 = 1 

222 PR. S2/10 

225 PR. ’Y', Y1 

230 FOR L = 1 TO S2 

240 Y6 = COS (Y1 * B(2,0,L)) 

250 Y7 = SIN (Y1 ^ B(2,0,L)) 

260 P4 = Y6 * G(L): P5 = Y7 V(l) 

270 D(L) = (P4+P5) 

275 T(L) = B(2,0,L) 

280 NEXT L 

n 

285 PI =: S2 

290 GO SUB 370 

300 X9 = (SQR (2/Pl)) * S5 

305 B( M, 0, Y2 ) = X9/1O00OO I 

310 X8 = EXP (-Y1) 
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312 B(7,0,Y2) = X8 

315 PR. X8, X9/I00000 

320 Y1 = Yl+0.1: Y2 « Y2+1 

330 IF Y1 < = ^01 GO TO 225 

340 S2 = 82 + 20 

350 IF S2t>=;^^0 TO 363 

355 M = Ivlfl ; PR, M 

360 GO. TO 2^ I Co 

363 STOP 

370 REIVi 

380 RBI MUAERICAL/INTEGRATIOM RCUTII4E 

390 REM 

400 E5=0: 05=0 

405 PR. »P1‘ , PI 

410 J5 = 2: P5 = P1--1 :N=1 

420 D8 = T(J5+1) -T(J5) 

430 F5 = D(l) + D(P1) 

440 IF N=(INT(N/2))*2 GO TO 470 
450 05 = 05+D(J5) 

460 GO TO 480 
470 E5=E5fD(J5) 

480 J5=J54-1 : N=N+1 

490 IF J5<=P5 GO TO 440 ; 

500 S5 = (F5f2*E5+4*05)*(D8/3,0) 

510 PR. ‘ INT’ , S5 
520 RETURN 
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This briefly explains the various quantities, which 
are shown in each of the plots given in the various chapters, 
A plot can display a single waveform as shown in Fig. 5.6(a) 
or (b). Or display two waveforms one superimposed upon the 
other as in Fig, 5, 2(a) or (b) or two waveforms, in which 
there is one displayed above and other below as depicted in 
Fig, 5, 8(a) or (b). In the case where only a single w'aveform 
is shown, the values are displayed above the plot and where 
more two waveforms are used, the quantities pertaining to 
both are separately shown, one above the plot and other below 
the plot. The values of these quantities are shown with the; 
help of a cursor and in the plots, the this is the third line 
For illustration, a typical plot in Fig. A is given, 
showing two waveforms, one above and other below. The plot 
atthe right, in particular, is chosen and the description f( 
each of the quantities is given. The six quantities above 
this plot refer to the upper signal waveform and three quant 
ties boiow the plot show the values for lower waveform. The 
description of these is as follows: 

a : This refers to the analysis operation done to get the 

waveform. In our case, since averaging was done to gej 

■ ' ] 

cleaner singal, acquired in channel l, therefore h 

’CHI AVERAGE* is printed in that place. 
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+HO rnimbor of such averagmgs done to 
b : This gives the numo^i oj. 

get the two waveforms. 

e : This xofors to tho arrpUtudo of ths waveform at the potot 
Whore cursor is placed. This is soon to bo 391. 

or mV. 

< the time instant on the waveform at which 
d : This refers to the tim 

" ’ the cursor refers to the value, in c. 

* ..ttho time instant of the plot and 
e ! This gives the starting 

always shows it as 000.00 msec. 

th total time to which the signal acquired rs 
f : This gives the tot 1 _ 

.. Thus > 2.4975 ms' means that signor 
displayed. Thus 2. ^ . 

only up to that time instant. 

r ■ s to the lower second signal, which was 

‘ ■ ”• '" 1 1 , 

acquired in channel z 

tho name ' CH2 AVERAGE' . 

ru omnlltude of the l<wei waveform at tho 
h ; This gives the ampi 

cursor point. . -o 

instant at which the value of h is given i 
i : The time ins tan v 

displayed in this. 

our description of the plots. 

This completes our 
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